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This paper aims to justify, at a microscopic level, the existence of a two-dimensional Bose metal, 
i.e. a metallic phase made out of Cooper pairs at T = 0. To this end, we consider the physics of 
quantum phase fluctuations in (granular) superconductors in the absence of disorder and emphasise 
the role of two order parameters in the problem, viz. phase order and charge order. We focus on the 
2-d Bose Hubbard model in the limit of very large fillings, i.e. a 2-d array of Josephson junctions. 
We find that the algebra of phase fluctuations is that of the Euclidean group E2 in this limit, and 
show that the model is equivalent to two coupled XY models in (2-|-l)-d, one corresponding to the 
phase degrees of freedom, and the other the charge degrees of freedom. The Bose metal, then, is the 
phase in which both these degrees of freedom are disordered(as a result of quantum frustration). We 
analyse the model in terms of its topological excitations and suggest that there is a strong indication 
that this state represents a surface of critical points, akin to the gapless spin liquid states. We find 
a remarkable consistency of this scenario with certain low-Tc thin film experiments. 
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I. INTRODUCTION 

The superconductor-insulator(SI) transition in low-Tc 
thin film systems [yj has drawn a lot of attention over 
the past couple of decades. These systems undergo tran- 
sition from superconductor (SC) to insulator as the dis- 
order, thickness or magnetic field is tuned. The problem 
has received a strong impetus after the experiment by 
Goldman et. al. |2] on homogeneous lead and bismuth 
films, which went from a superconducting phase to an 
insulating phase as a function of thickness, and which, 
right at the interface, was probably metallic. Since elec- 
trons usually do not form a metallic state in 2-d, it was 
argued that [Ref.Sb] this T — transition is due to the lo- 
calisation of preformed Cooper pairs. It was also claimed 
that the resistivity at the transition is univcrsal[Rcf.3b]. 
Similar SI transitions have been observed in granular su- 
perconductors and Josephson Junction arrays [El. This 
scenario has been called into question after a recent mag- 
netic field tuned experiment in the Mo-Ge sample [H, 
where the metal is no more a point in the phase diagram, 
but exists as a separate phase. We would like to point 
out that this is not the first observation of a metallic 
phase in a two-dimensional, otherwise superconducting, 
system. We came across, at least two separate instances 
of this phenomenon in granular superconductors — one 
in Ga films H and the other in granular Pb films P, 
where a metallic phase is found to be sandwiched be- 
tween the superconducting and the insulating phases. A 
similar observation has been reported in Josephson Junc- 
tion arrays M . Each of these observations has probably 
a detailed explanation within the scope of the specific 
experimental system being measured. However, there is 
something common in these systems which is quite hard 
to overlook, viz. the preformed Cooper pairs are very 
much alive when the metallic phase is seen; and continue 



to exist till the insulating transition and beyond. Moti- 
vated by this fact, we would like to ask the question — is 
it possible for the charged bosons, i.e. Cooper pairs, to 
form an incoherent metallic phase at T = 0? The intrigu- 
ing feature of this phase is that although the bosons are 
mobile, they do not bose condense at any temperature. 
Instead, they are dissipative and fail to drive a supercur- 
rent even at T = 0, unlike a superconductor. We regard 
this phase as a Bose metal(BM). In this paper, we give 
arguments justifying the existence of a Bose metal in a 
physically realizable system. 

The natural question, then is, what has been missing 
in the current theoretical models where a metallic phase 
was not obtained. Our thought on this issue is the fol- 
lowing: most of the earlier theories [pPfl tried to attack 
this problem from the superconducting side of the phase 
diagram and projected it onto a basis diagonal in the 
phase states. In all these theories, there was a single or- 
der parameter, viz. ip —< e**^ >, where (j) refers to the 
phase of the charge boson. In the superconducting phase, 
this order parameter is well developed. At the SI tran- 
sition, tp ^ and the phases of the charge bosons get 
scrambled. The scrambling of the phases y,|l2l had so 
far been taken as the indication of onset of the insulating 
phase. However, this is not enough to characterise the 
insulator. A Bose insulator (BI) phase is characterised by 
an extra order parameter, viz. the charge density. It is 
like a charge density wave but built out of Cooper pairs. 
This piece of physics has been missing in the existing the- 
ories P,p|-[lO[ . The central point of this paper is that the 
phase fluctuation physics of superconductors should be 
viewed as a two order parameter problem, viz. there is 
a competition between phase order and charge order. It 
almost follows from this fact that the destruction of one 
order parameter does not necessitate the growth of the 
other order parameter. This implies a possible existence 



of a disordered phase where both the order parameters 
are zero at T = 0. We consider this to be a Bose metal 
phase. 

Secondly, we work in the limit when the average filling 
of the bosons(rio) per site is very large. In this limit, 
the Bose Hubbard model [Ref.3b], which plays a central 
role in the SI transition, becomes equivalent to a Joseph- 
son junction array (JJA) model. Most of the existing 
work on the JJA model does not treat this large no limit 
consistently. In many cases, the constraint due to the 
average filling, eqn.(2) below is neglected [Q; in certain 
other cases, there is a tendency to replace the E2 alge- 
bra, which is the appropriate algebra for this limit (see 
sections II and V for more details), by a qualitatively 
different algebra, viz. SU{2) algebra, which holds in the 
hard core limit, i.e. for small fillings. Hence, we shall 
give a quite detailed description of the basic formulation 
of this limit in this paper. 

Thus, we consider a pure model of Josephson-coupled 
Cooper pairs with extremely large fillings, interacting via 
on-site and near neighbour repulsive forces in two dimen- 
sions. This model captures the basic physics of granular 
superconductors Q, except that disorder is absent. It 
is applicable to a limited extent to the Josephson junc- 
tion arrays as well. We are able to demonstrate that the 
physics of phase fluctuations, within such a model, can 
be described by two coupled anisotropic XY models in 
(2-f l)-d, one corresponding to the charge degrees of free- 
dom and the other the phase degrees of freedom, where 
the coupling is "XY-like" . The two mechanisms which 
drive the two transitions are as follows - disordering of 
the vortices and the bose condensation of vortices. Phase 
order is destroyed when the vortex-antivortex pairs un- 
bind(or, vortex loops blow up in (2-|-l)-d) as a result of 
the quantum fluctuations; the charge order grows when 
the vortices Bose condense [|ll[ . In light of this, a search 
for a completely disordered phase in the charge picture 
translates into a search for a non superfluid(SF) liquid 
in the vortex picture. The presence of unbound, as yet 
uncondensed, dissipative vortices makes the disordered 
phase metallic. The two mechanisms mentioned above 
are separate processes owing to the presence of retar- 
dation, or equivalently, dissipative effects, coming from 
a gauge field mediated interaction in the vortex picture 
(for more details see sections lid and Illb). In the charge 
picture, the existence of a disordered phase like that of 
a Bose metal results from quantum frustration effects, 
i.e. the zero point motion is not strong enough to set 
up superconducting correlations and the long range in- 
teractions are not sufficient to set up a charge ordered 
state. This state seems to be intrinsically related to the 
(gapless) spin liquid states. Although disorder is not ex- 
plicitly present in our model, our results on the nature 
of the Bose metal phase can be readily generalised to 
the case where disorder is present. We find a remark- 
able agreement of our predictions with the Gallium film 
experiment |p[ . We believe this provides evidence of the 
existence of such a phase and supports the scenario de- 



scribed above. All of our considerations in this paper 
are restricted to T = and zero magnetic field unless 
mentioned otherwise. This situation is relevant to the 
experiments in Refs. Q and g. 

Before we go over to the main part of the paper, we 
would like to mention that the idea of two order pa- 
rameters is not completely new and has appeared in the 
discussion of the hard-core limit of the Bose-Hubbard 
model [E2| and also in the gauge theory description of 
the Josephson Junction arrays ||l3| . The problem with 
the former is that superconductivity is never destroyed 
away from half-filling in their model (please see section 
lib), while the latter invokes the idea of self-duality which 
implies logarithmic interaction among Cooper pairs and 
is not a good representation of the realistic thin film sam- 
ples. 

Apart from the existence of a novel disordered Bose 
liquid state, there is another important feature which 
emerges from our work. It is well known that |14[| non- 
interacting electrons cannot support a metallic state at 
low temperatures in 2-D in the absence of spin-orbit scat- 
tering. The situation is much less clear when both in- 
teractions and disorder are present. This fact has been 
brought into the limelight with the recent observation 11^] 
of metal-like states in two dimensional electron and hole 
systems in semiconductor based materials. In the wake 
of these observations, our investigations suggest that a 
metallic state is a possibility if the electrons bind them- 
selves into Cooper pairs and behave as if they arc bosons. 

The plan of the paper is as follows. In section II, we 
introduce the model and straighten out some of the basic 
issues which relate to the model. Using a combination of 
qualitative and quantitative arguments, we explain why 
the BM phase is feasible. The mapping onto the cou- 
pled XY models is demonstrated here. We quantify the 
arguments in section III and sketch the phase diagram 
of the model. Section IV contains an estimate of the re- 
sistivity of the Bose metal phase and a comparison with 
the Gallium film experiment. The thermodynamics of 
this strange metallic state and its connection with the 
spin liquids are discussed in section V. We wind up with 
a discussion of certain relevant issues and conclusions in 
section VI. Some of the calculational details may be found 
in the appendices. 



II. THE FOUNDATIONS 

A. The model 

We shall consider the following model Hamiltonian in 
this paper — 

iJ = - J ^ cos{(j)i - (f>i+a) + Vb ^((5rii)^ 

+ Vi'^ {Sn.i + 5ni+af - p-^Sui- fj,Nno (1) 



where drii = rii — no, with rii = number density operator, 
no = neutrahsing background charge density(or, equiva- 
lently, the average density of Cooper pairs), N = number 
of lattice sites, and jl — ^-^Vo the renormahsed chemical 
potential in the problem. Also, a — x,y refers to the spa- 
tial unit vectors. The first term in eqn.(l) represents the 
kinetic energy of the bosons(Coopcr pairs for our case), 
the second term the onsite repulsion among them (which 
should be non-zero to prevent any collapse of the bosons 
onto a single site) , the next term the repulsion among the 
nearest neighbours which acts to set up a charge order 
in the system, and the last two terms are the chemical 
potential terms. Eqn.(l) has to be supplemented with 
the constraint -^ Si ^ '^^ ^~ "-Oi which in terms of the 
charge fluctuation operators convert into 



[p,pt]~o 



(6c) 



y^ < Sui 



>=0 



(2) 



This equation implies /l = and, in the rest of the pa- 
per, we shall forget about the chemical potential terms 
completely. 

On the experimental side, the three parameters of the 
Hamiltonian (1) can be determined in the following way 



J = (i?Q/2i?„)Ao 



(3) 



Here Rq = [h/Ae^) = QAhKQ., Rn = normal resistance 
of the film, and Ao = pairing gap. The interaction con- 
stants Vq and Vi are related to the inverse of the capaci- 
tance matrix C^ of the grains. 



B. The commutation relations 

The model discussed in eqn(l) needs to be supple- 
mented with the phase fluctuation algebra, which con- 
stitutes the appropriate commutation relations for this 
problem: 



[5ni,<j)j] = iSi^ 



(4) 



Eqn.(4) implies an angular momentum representa- 
tion[Ref.3a] 



5n,, = 



d 

d(j)i' 



Further, let us define the operators — 

L = i — , andP^e^-^. 
d(p 

Thus, eqn.(4) can be recast as — 

[L,P] = -P 



(5) 



(6a) 



with PP' = I. It is well known in quantum optics that 
the phase operators being ladder operators [as is seen 
from eqns.(6a) and (6b)] usually do not commute p7|. 
However, in the large no— limit, they do (see appendix 
F). Thus, eqn.(6c) is strictly valid in the large no— limit 
of the problem and, hence, our discussion holds good in 
this limit only. Now, eqns.(6) constitute the algebra of 
the Euclidean group E2 [Q , the group of translations and 
rotations in 2-d, with the square of linear momentum(P) 
restricted to unity for our case. It deserves to be men- 
tioned here that the SU{2) algebra used in the context of 
the hard-core limit of the Bose problem |l2| is distinctly 
different from this algebra and is obtained in the opposite 
limit, i.e. the small no— limit of the Bose Hubbard model. 
The effects of the change in the group structure are quite 
significant. In the hardcore model, away from half fill- 
ing, increased interactions changes the SF to a supersolid 
phase. Thus, superfluidity is never destroyed away from 
half filling. On the other hand, our model allows super- 
fluidity to be quenched at arbitrary fillings. Further, the 
conserved quantities supported by the two algebras are 
different: the invariant of E2 algebra is (square of linear 
momentum) P^ + Py = p^ (= 1 for our case), as com- 
pared to Sl + S'^ + S'^ = s{s + 1) of the SU{2) algebra. 
Thus, unlike the SU{2) case, (a) here the constraint on 
the z component of the spin is much weaker and (b) the 
z component of the spin enters very anisotropically in 
the algebra compared to the x and y components of the 
spin. (Here operator L is referred to as the z component 
of spin for E2 algebra; see section V for more details.) 
Hence, in the disordered/symmetric phase, although the 
X and y components of E2 spin might acquire a gap like 
the SU{2) spin, the z component of the former might re- 
main gapless unlike the latter. Put in simple words, there 
may be a length scale determining the local superfluid- 
ity in this phase and as yet no length scale associated 
with the charge ordering in the system. A calculation 
of the charge charge correlation function indeed justifies 
this, as discussed in section V. Thus, given a disordered 
phase within a model built out of E2 spin operators, it 
can support gapless excitations. We shall give detailed 
arguments in sections IID and HI as to why the model (1) 
along with commutation rules (6) contain a completely 
disordered phase. 

Now, if we rotate the charge fiuctuations Li at each 
site, viz. 

Li -^ e"5-''*Li, where Q — (7r,7r) 

we obtain the Hamiltonian as 

H^Y. ■^'^PkPk + J2 VkLlLk, (7) 

A; k 

where 



[L,pt]=pt 



(6b) 



Jfe = —J{coskx + cosky) 



Vk = (Fo + Wi) - 2Vi(cosk^ + cosky). 

and P^L couple through the commutation relations (6). 
In this form, it assumes the shape of a 2 order parameter 
problem. (Macroscopic occupation of A: = mode of P 
reflects superconductivity and that of L a charge density 
wave.) However, because of commutation rules (6), it is 
very hard to diagonalise this hamiltonian in this form. 
So, we seek alternative means. 



C. The charge picture 

In this section, we demonstrate that the model Hamil- 
tonian (1) is equivalent to two coupled XY models in 
(2+l)-d. To do this, we first write the model (1) in terms 
of a path integral representation |jlq,tt9[ keeping the com- 
mutation relations (6) in mind. The partition function 
Z = Tre~^ , where /3 is the inverse temperature {h and 
ks are taken to be one, unless mentioned otherwise), can 
be written as — 



(l)We decouple the (V^mi)^ term using a Hubbard- 
Stratanovich field pia, viz. 



/ PK„e-(^/^^^^^) S.c ^'"-^ S. ™'(v. 



P.c) 



(2)First, we notice that the coupling term with rrii is 
invariant under shifts of pia by 2tt. So, we break up the 
intgral over pta from — oo to -l-oo into that of periods of 
27r. Further, we split up pia into a curl and a gradient 
part. Since the divergence oi pia couples to rrii, only the 
gradient part enters the dynamics. Thus, we obtain for 
this part. 



(l/2yiAT)^. cos(V„0i) -i^iTiiV^ 



Z = 



E 

{mi{T)} 



27r 



V(j),{T)e- 



where 



V0,e 



where we have used an (inverse) Villain transformation 
p3| . (3) Now, using an (inverse) Villain transformation 
again, one can execute the sum over integers rui — 

^ exp(-i^m,(e-'Q-'-'V,0, -K V^e,) - FoAr^m^) 

{rrii} i i 



+ Fo E ™?(^) + ^1 E('"'(^) + ™«+"(^))'] (8) 



Putting all these together, one obtains 



27r 



v<t>,ve,e- 



with 



with periodic boundary conditions in the imaginary time 
direction implied. Here a — x,y refers to the nearest 
neighbours in the space direction and miij) are integers 
meaning the change in the number of Cooper pairs from 
the average at the site i. We shall consider only the 
case T = 0, so that the integral over imaginary time 
extends to infinity. Next, we discretize the time axis 
with an interval Ar and rescale the lattice constant in 
the space directions to be unity. Further, we rotate the 



5 — J At y^ cosCVaCJ^i) 



'Wr(f>i+^ 0, 



(10) 



integers nii^T) 
we obtain — 



-iQ-Ti 



mi{T), where Q = (tt, tt). Thus, 



i i 

-t- ViAt^(VqTOj)^ - JAr^cos(VQ(/>j). (9) 

ia ia 

Here summation over i refers to the time axis as well, 
and so does the index i in rrii and 0^. The derivatives 
V^ in eqn.(9) and also in what follows are lattice deriva- 
tives [po| . To avoid any confusion, we shall reserve the 
notation r^ for the spatial coordinates of the i*'* point 
and Xi its space-time coordinates. 

To show that the action (9) is equivalent to 2 coupled 
XY models, we follow the following sequence of steps. 



where the phase (f> is associated with superfluidity and 
the phase 9 with charge density wave. Eqn.(lO) explic- 
itly shows that the phase fluctuation physics of (granular) 
superconductors, in the absence of disorder, is equivalent 
to two coupled XY models in (2-|-l)-d, where the cou- 
pling is XY-like. A few comments are in order. In the 
limit Vi = 0, we have 9i = 0, the 6'-terms in action (10) 
drop out, and we obtain a single XY model in (2+1)- 
d, as has been discussed previously P, |l0| , |l9| . Also, the 
coupling term (the second term) is highly anistropic in 
(/) and 6, explicitly breaking self-duality in this system, 
assumed in Ref. ]13]. Eqn.(lO) is one of the key results 
of this paper. Thus, we see that the destruction of su- 
perfluid state is driven by one XY model, whereas the 
other XY model characterises the growth of the charge 
ordered state. And, hence, as the parameters are tuned 
one transition does not necessarily accompany the other. 



D. The vortex picture 

Although the action (10) shows that the superfluid or- 
der and charge order are driven by different XY mod- 
els, it is not clear whether the coupling between them, 
which is quite complicated, guarantees a completely dis- 
ordered phase. In order to answer such a question, we 
now consider the model (1) in the dual picture, i.e. of the 
vortices. To do this, we invoke a duality transformation 
P,OyOJl9]. Starting from action(9) pil, one can show 
that 






l<7h) 



T ' q,uj 






(11) 



where j^ ^^ are Fourier transforms of the integer vortex 
variables jf (j" — vortex density, jf — vortex current), 
and 



G-\u,q)^[u;^ + 






FiAr 



ln{^ 



■fiq-Qf 



J At' 



(12) 



The intervening steps are quite standard 123] and are 
discussed in Appendix A. w, q refer to lattice frequency 
and momentum respectively [Appendix A]. One may like 
to note here that there is no Magnus force term on the 
vortices, corresponding to the disappearance of riQ from 
the problem in the large uq limit. 

In the previous subsection, we observed that when 
Vi — > 0, the model is equivalent to a single XY model 
rather than two coupled XY models. The question, obvi- 
ously is, how this change is reflected in the dual picture. 
The basic answer lies in the appearance or disappearance 
of the zero point motion term for the vortices (more pop- 
ularly known as the vortex mass term). To show this, 
we look at the long wavelength low frequency modes. In 
the limit Lu,q —> and to <C c^g (where Cg is the plasmon 
velocity), the Green's function G{uj,q) splits up into a 
singular part Gs{uj,q) and a constant part Gq, viz. 



G{LU,q) ~ Gs{uJ,q) + Go, 



(13) 



where 



Gs{iO, q) = 1/(C^2 ^ ^2^2^^ g^ ^ ^2/8(^2 ^ ^2)2 



with 






VqAt 



9= ""^ V fO —jO 

;„( 2 )Z^J«,.. 2-^-9.— 
yjAr^ q.u] ^ 

+ TT-^^T.^U^si^^'i) + Go]j%.-^, (14) 

Since Go is a constant, it allows us to identify this term as 
the vortex kinetic energy, i.e. mass term. From (14) and 
(13), we notice that this term does not exist when Vi — 0. 
Thus, we have shown that the change of the nature of the 
XY models is tied to the existence or non-existence of a 
vortex mass term. 

We now look at the other terms in the action (14). 
The first term in the action is the usual logarithmic in- 
teraction term among the vortices, and the second term, 
because of the retardation effects, leads to a dissipative 
term for the vortices [|l6[. In the limit | r.i — tj \<^ Cs \ 
T ~ t' \ (which is the same quasistatic limit discussed so 
far) , this part of the action takes the form — 



c2;n(2/JAr) 



i,3 



Y^LoHni- 



1 



where qi = ±1 refers to the charge on the i vortex. 
Physically, the aforesaid limit corresponds to the slow 
motion of vortices. The source of this heat bath(or, dis- 
sipation) is the gauge field (discussed in Appendix A), 
or more precisely, the transverse modes arising from the 
quantum fluctuations in the system. (There will be addi- 
tional retardation effects in a real system from external 
heat bath mechanisms, e.g. electrons in the vortex core, 
etc.) These features of vortices in granular supercon- 
ductors have been discussed previously by Eckern and 
Schmid !§. 

One may also undo some of the steps in Appendix A, 
and recast eqn.(14) in terms of a gauge field with appro- 
priate action as — 



2^^Y,J^A^+^^^ 



Gn 



Hihf) 



/ J J q,ujJ — 



q^ — LU 



(15) 



Thus the action is 



(with the gauge condition Vq.^" = imposed). A few 
comments are in order. Eqn. (15) shows that, considered 
in the vortex picture, the quantum phase fluctuations 
in a 2-d superconductor, as described by model (1), is 
equivalent to a two-component quantum plasma (bosons 
of two flavours, viz. vortices and antivortices) moving in 
a fluctuating gauge field A^^ Q . Secondly, the above sce- 
nario is a simple quantum mechanical extension of clas- 
sical phase fluctuations in a 2-d superconductor, which 
is described by a two component classical plasma under- 
going screening by a static electric field [E = — V^o), 
as described by Kosterlitz and Thouless(KT) IHl. The 
effect of including quantum mechanics in the problem, 
apart from bringing up the importance of the quantum 



statistics of vortices, is to make the electric field dynam- 
ical, viz. E = —VAq — {l/cs)dA/dT with a magnetic 
field B — W X A. Whereas the importance of the statis- 
tics is to allow for the superfluidity of the vortices, an 
important consequence of the dynamical nature of the 
electromagnetic field is that there are retardation effects, 
viz. 



E 



drdr'- 



n{T).rj{T') 



^{n{T)-r,{T')r+cl{T-T') 



2' 



which break Galilean invariance. (This feature is inten- 
sified in a real system by external heat bath mechanisms 
mentioned before pq| .) This is not surprising because 
the action (15) has the structure of Maxwell's action, 
which is reputed to have Lorentz invariance but lacks 
Galilean invariance. Now the absence of Galilean invari- 
ance will have a strong effect on our system, because 
it is bosonic. It is well known that all the delocalised 
bosons will condense into the supcrfluid state only if the 
system is Galilean invariant ||26| . If this invariance is ab- 
sent in a Bose system, then as the relevant parameter is 
tuned, the (vortex) condensate is gradually depleted and 
at one point the superfluidity will be completely lost. 
At this stage, it is important to recollect what impor- 
tant processes are going on in this system. There are 
two of them: destruction of vortex (and antivortex) su- 
perfluidity (owing to retardation effects) and the binding 
of vortex-antivortex pairs, corresponding respectively to 
the destruction of charge order and the growth of phase 
order. Now, these two processes are controlled effectively 
by two separate parameters, viz. c^/g^ ~ (Vq -I- 8Vi)At 
(corresponding to the strength of the retardation effects) 
and g"^ ~ J At (corresponding to the strength of the log- 
arithmic interaction) respectively, where At ~ 1/Ao pj] 
(please refer to appendix D and section III for the no- 
tation and appropriate details). As a result, the vor- 
tices(and antivortices) do not necessarily condense into 
the supcrfluid state as soon as they unbind. This leads 
to the possibility of a non-SF vortex liquid, or equiva- 
lently, a BM phase. In fact, that's what we find when 
we quantitatively evaluate these processes in section III. 
There is a simpler way to see what is happening here. As 
we noted earlier, following eqn.(14), the kinetic energy of 
vortices originates from two sources — (a) quantum zero 
point motion (the Gq term) and (b)action of an effective 
heat bath (the Gs term). In the delocalised state, if the 
source (a) dominates, the vortices move coherently, and 
since they are bosons, they form a supcrfluid. On the 
other hand, when the source (b) dominates, because of 
the random nature of the effective heat bath, the motion 
of the vortices is necessarily incoherent, and one is in a 
metallic phase. In this phase, no charge order is set up 
and one obtains a Bose metal phase. 

These features of dual vortices are not special to a lat- 
tice model, but observable in the continuum formulation 
as well, as discussed in Appendix B. 



III. THE PHASE DIAGRAM 

In the previous section, we argued why the Hamil- 
tonian described by eqn.(l) may contain an incoherent 
metallic phase. In this section, we quantify these argu- 
ments by calculating the phase diagram of model (1): we 
shall locate the phase boundary where superconductiv- 
ity is destroyed and the one where charge order is estab- 
lished. 



A. Destruction of charge superfluidity 

This is done non-perturbatively by estimating where 
the vortex loops blow up in (2-|-l)-d [g7|. This happens 
when the entropy of the loops overcomes their interaction 
energy J28| . A good estimate of the interaction energy is 
the self energy of the loops, simply because dipole-dipole 
interactions fall off as l/r^ and the mutual interaction 
energy of the links in a loop is much smaller than the 
self-energy of the loops when the loops are fairly large. 
Thus, the effective free energy of the loops is given by 



J^ = 



Hjh) 



G{0)-fii]N 



(16) 



where N is the number of the links in a vortex 
loop, G(0) — diagonal part of the Green's function 

~ /- 1^ i'2rf ^(^''^)' obtained from cqn.(12), and /i; = 
entropy of the loops = ln3 for our case [g9[ . An estimate 
of G(0) is given in the Appendix C. The loops blow up 
when J^ < 0. Thus, superconductivity is destroyed when 



(Vo/J) + 8{Vi/J) > 6o, 



(17) 



where bo — 2 (69//^;)^ ^^^ bQ = a number of order unity ~ 
3.1725, defined by eqn. (C6) in appendix C. This phase 
boundary has the character of (2+l)-d XY model, at 
least when Vi <^ Vq (as seen from eqn. (10)), and hence, 
the superfluid density changes continuously across this 
transition line. 



B. Destruction of vortex superfluidity 

It is well established that the existence of charge order 
implies superfluidity of vortices and vice versa p^ . We 
shall follow this notion here and estimate the growth of 
charge order in terms of the superfluidity of vortices. We 
mentioned in section Ild that the vortices and antivor- 
tices move in the presence of a dynamical gauge field and 
argued that because of the latter there are retardation ef- 
fects which deplete the (vortex) Bose condensate (due to 
broken Galilean invariance). As a result, the vortices and 
antivortices do not necessarily condense into a superfluid 



state as soon as they unbind. This physics of suppres- 
sion of Bose condensation as a result of gauge field fluc- 
tuations is not new but has been explored substantially 
in the context of spin charge separation theories in high 
Tc superconductors |31^,G3]. The discussion here is very 
similar in spirit to that piece of work. 

To estimate the strength of the parameters where 
vortex superfluidity is destroyed, we follow the self- 
consistent functional approach of loffe et. al. psl. They 
did the calculation for a very similar piece of physics, i.e. 
how the coupling to a gauge field can kill supcrfiuidity 
in a bosonic system (with logarithmic interaction) as a 
result of broken Galilean invariance stemming from retar- 
dation effects. We refer the reader to that paper for a full 
description of this technique. A short discussion of this 
is given in Appendix D. We can follow their approach 
here, simply because the lattice action and the contin- 
uum action have identical structure in the low frequency 
long wavelength limit, as discussed in Appendix B. To do 
this, we first replace the two component plasma by a one 
component plasma, i.e. charges(vortices) of one flavour 
moving in the background of flxed neutralising charges of 
the other flavour. This approximation is very standard 
and captures the salient features of the problem, until 
and unless the plasma is extremely dense 32 1. Then, 
the results of sec.V in Ref. |2^ can be directly carried 
over here with the identification of the Coulomb inter- 
action parameter ac and transverse gauge field coupling 
constant ag 
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ac = (^V2)«i/(wo + 8wi)2 
ag = (l/7m„)[H-(uo/8ui)] 



(18) 



where n„ = average vortex density and Vi = {Vi/J),{i ~ 
0,1). Since what counts in the destruction of vortex 
superfluidity are the free vortices(and antivortices) , we 
take riy ~ nf, which goes inversely as the square of 
the correlation length f+ which diverges at the SC- 
BM boundary from the BM side. Thus, from (17), 
rif ^ [feoC^^o + &Vi)/J - 1]2^ {bo = l/bo) with u ~ 2/3, 
since this phase boundary has the character of (2-|-l)-d 
XY model. Also, if we are above the phase boundary of 
eqn.(17) (see fig.l), we have a <C 1 as in Ref. |g^, where 
a = yJoTc- Since as we shall note below that the transi- 
tion from SF to non-SF state takes place at a small value 
of a, this assumption of small a is self-consistent. 

Let us first consider the simple case when we are 
slightly above the phase boundary (16) and along the 
V\ axis (see fig.l), so that Vb = 0, Vi ^ 2.1J. Then we 
have, a ~ 0.2 and a^ ^ 1; the calculation of Ref. |2^, in 
that case, suggests we are in the disordered BM phase. 
Thus, there is at least a small region close to the phase 
boundary (17), where the system is metaUic. 

We now complete the calculation of phase boundary of 
BM-BI transition, following Ref. Q. Since < rii, < 1, 
from (18), (Xg > 0.32. From fig.7 of Ref. |||, one can 
see that in this regime, the phase boundary tends to sat- 
urate at a = ttcr — 0.08. Actually, for larger values of 



ag, acr is a little less; but, the point is that acr is always 
finite, however small it be, for large values of ag. That is, 
the vortices always form a superfluid for small enough a. 
A physical way of seeing why there is always superfluid- 
ity for small enough Coulomb repulsion is as follows. The 
strength of the Coulomb repulsion is controlled by the ef- 
fective charge g (see appendix D), which also controls the 
strength of the coupling between the particles (vortices) 
and the transverse gauge fleld. So, when the magnitude 
of Coulomb repulsion is small, g is small, and hence, the 
coupling with the transverse gauge field is small as well. 
As a result, the particles do not feel the effect of retar- 
dation strongly enough in this limit and condense into a 
superfluid phase. Thus, the vortices are in a superfluid 
state when a < acr, which means — 



1 



(SVi/J) 



[{Vo/J) + {8Vi/J)f > Co 



(19) 



where cq = 4(7r/8Qcr)^ — 96 for a^ — 0.08. This calcu- 
lation implies a jump in the (vortex) superfluid density 
at this phase boundary [ p3[ . This means that the phase 
transition is either of flrst order or has a KT character. 
More calculations are necessary to resolve this point. 

We display a schematic phase diagram determined by 
the eqns.(17) and (19) in flg.l. Eqn.(17) is represented 
by the curve LXM and (19) by ZYX. They seem to meet 
at a tricritical point X. Since the retardation effects are 
very strong in this model (meaning that the value of acr 
is small), the constant cq in eqn.(19) is roughly an order 
of magnitude larger than the constant feg appearing in 
eqn.(17); and hence we expect that this crossing between 
the two curves will always occur, implying the existence 
of the Bose metal (see section VI as well). The vari- 
ous phases determined by these eqns. are as marked in 
fig.l — when the zero point motion of the Cooper pairs 
is large, the system is superconducting (region LXMO); 
when the interactions dominate, the system is charge or- 
dered and insulating as a result — we call it a Bose in- 
sulator (BI); and in the intermediate region, the system 
is disordered and, hence, metallic (please see section IV 
below) — the BM phase (region ZYXL). 

What one sees on the phase diagram of fig.l is that 
the metallic phase is more prominent towards the Vi axis 
rather than the Vq axis. Further, from eqn.(l) we note 
that the Vi-term contains both on-site and nearest neigh- 
bour repulsion energies. This shows that the Bose metal- 
lic phase is to be expected in cases where these energy 
scales are of comparable order of magnitude, a situation 
well represented by the granular superconductors 0. 

We can provide a physical explanation in the charge 
picture as to why the metallic phase opens up along the 
Vi axis: let us focus on eqn.(l) and say that supercon- 
ductivity is already quenched. First consider the case 
Vo large and Vi small, so that we are along (or, close 
to) the Vo axis. Naively, one would expect from the Vi 
term that Srii + Sui^a ~ large. But, this costs a large 
energy from the Vq term. As a result, what is favoured 



is Sui ~ 0, i.e. n^ ~ no, a (rather trivial) charge or- 
dered state. Now, consider the opposite case: Vi fairly 
large and Vb small, so that we are close to (or, along) 
the Vi-axis. Since Vi is appreciable, eqn.(l) suggests 
Srii + Sui^a — 0. A non-trivial configuration may be 
Srii = —SriiJ^a = 1, for some z's and zero otherwise. This 
state represents an RVB-like state of fluctuating charge 
(particle-hole) dipoles, the equivalent of spin singlets here 
(please see section V for more on this point). Vq being 
small, the VoJ2ii^''^i)'^ term does not cost much energy 
for this kind of state. This is the disordered BM phase. 
When Vi is increased further, the dipoles freeze into a 
charge ordered solid. Thus, very crudely speaking, the 
smallness of Vb is a source of frustration in this model for 
finite Vi. 

Since J is inversely proportional to the normal resis- 
tance Rn [eqn.(3)], as the normal resistance of a thin film 
is tuned in an experiment we gradually cross from super- 
conducting to non-superconducting regions according to 
the phase diagram of fig. 1. Some such typical traces are 
shown in the figure — the dashed lines OA, OXB and 
OC. Traces OXB and OC represent a superconductor - 
insulator transition, a case which has been discussed very 
widely in the literature 0] . On the other hand, the trace 
OA represents a superconductor - metal - insulator tran- 
sition. In this part of the phase diagram, the retardation 
effects are very strong and the system passes through an 
intermediate disordered phase. This is a new prediction 
made by our analysis. Model (1) which leads to this is 
good for s-i-s junctions. 



IV. RESISTIVITY OF THE METALLIC PHASE 

The discussion in the previous section suggests that 
the appropriate model for the Bose metal phase is that 
of uncondensed bosons (vortices and antivortices) in a 
transverse gauge field. In the frustrated BM phase, the 
vortices are unbound, but they fail to bose condense due 
to retardation, or equivalently dissipative, effects. As 
a result, the system exhibits metallic behaviour. This 
section is devoted to making a simplest possible estimate 
of the resistivity of this metallic phase. We shall evaluate 
the resistivity of the charge bosons pc in terms of the 
conductivity a^ of the dual variables, viz. the vortices. 
The relation between the two is given by (we mention 
the factor of h explicitly in this formula) — 



{h/4e^ 



(20) 



Also, ffy is given by the Drude formula ay = nyfTtr/my, 
where Ttr refers to the transport time and n„/ is the 
free vortex density. There are three contributions to the 
conductivity of the vortices — (i) dissipation due to the 
action of a heat bath(coming from the trasverse modes), 
(ii) scattering from impurities, and (iii) Bardcen-Stcphen 
processes. In this part, we shall focus on the contribu- 
tion from process (i). We give here a short description 



of how this dissipation mechanism comes about in this 
BM phase. As one enters this phase from the SC phase, 
say (following curve A in fig.l, for example), the vortices 
and antivortices unbind in 2-d. These unbound charges 
screen each other leading to a finite screening length, 
i.e. a gap in the longitudinal part of the gauge field. 
However, the transverse part of the gauge field is gap- 
less, because there is no spontaneous symmetry breaking. 
These transverse modes have their source in quantum 
fluctuations and mediate the dissipation process. They 
represent the plasmons. (Physically, this makes sense, 
since the plasmons are gapless in 2-d.) The physics is 
essentially that of damping which results when a charge 
(vortex) moves through a background of other charged 
particles. This lends a hydrodynamical character to this 
dissipation mechanism. In section lid, we saw that the 
most important modes are uj <ti Cgk. So, we shall con- 
sider contributions from these modes only. The leading 
order contribution, then, comes from the a; = 0, /c = fi- 
nite modes, and we make an estimate from this sector 
only. We think the higher order contributions will lead 
to simple renormalisation of coefficients, so long as the 
density of vortices (and antivortices) is low (please see 
below). Thus, we are led to evaluate the resistivity of 
a set of uncondensed bosons moving in a static random 
(but, annealed) gauge field, with a variance 

< H{q)H{-q) >= l/2c^Jn{2/JAT) 

as seen from eqn.(15), where H = W x A. Now, also in 
this delocalised phase, the interactions between the vor- 
tices are screened (by antivortices) and, in a lowest order 
approximation, we shall treat them as non-interacting 
bosons. Thus, the transport time is the same as that 
of a non-interacting particle scattering from a random 
magnetic field, which is given by ||33|| — 



Ttr 



■^.2{Vo + 8Vi)At 



where we have used the above field distribution and the 
parameters as enumerated following eqn. (13) . The lattice 
constant in the time direction. At ^ I/Aq |3J]. Thus, 
we have — 



;RQ{nyfeo){VQ + 8Vi)/Ao 



(21) 



where Rq = h/Ae^ = QAbKQ,. The quantity n^f^Q is 
inversely proportional to the square of the superconduct- 
ing coherence length which diverges at the SC-BM phase 
boundary, i.e. n./^g _ (^o/^^)2 _ [{Vo + Wi)/k,J~lf\ 
obtained from eqn. (17). Using eqn. (3), we can rewrite 
this as — 



nvf^l 



^Rc ' ' 



(22) 



where Re = boRg/KVo + 81/l)/Ao] is the critical resis- 
tance of the film where the metallic phase sets in and 



V = correlation length exponent. The exponent v is de- 
pendent on the particular universality class of SC-BM 
phase transition. For our case, the SC-BM phase bound- 
ary has the character of a (2-|-l)-d XY model and, hence, 

V ~ 2/3. However, a real material, like the low- Tc thin 
film systems, is heavily disordered and v will certainly be 
very different from this, as we shall see below. So, from 
eqn.(21), we have 



RuIRq ~ feo[(i?n/-Rc) - if^'iRqlRc) 



(23) 



The case when disorder is present — The foregoing dis- 
cussion can be readily generalised to the case when dis- 
order is present. Owing to the Drude formula, the struc- 
ture of eqn.(21) is not unique to model (1) but will come 
about when the bosons form a metallic phase as a result 
of phase fluctuations (in zero magnetic field), i.e. 



Rn 



2RQ{nyf£,l)fii 



(24) 



where ^d represents the vortex mobility, or the fric- 
tion factor, arising from the dissipation of vortices and 
is model dependent ||3§]. For the pure case discussed 
above, the damping is due to the transverse modes and 
liD ~ (Vb -|- 8Vi)/7r^Ao. When the normal resistance is 
tuned through R^ as we saw above for the pure case, 
the most important dependence of BM resistivity in the 
formula (24) comes from the free vortex density ni;/Co, 
i.e. 



^vf^O 



^Rc ' 



where the critical exponent v characterising the SC-BM 
phase transition depends on the particular universality 
class being considered. As compared to this, the dis- 
sipation factor /^£), in general being dependent on the 
morphological and normal properties of the film, is only 
weakly dependent on the normal resistivity about Re [p5[ . 
Thus, in general. 



Rn 



RQ[iRn/Rc)~l] 



2v 



(25) 



For the case of Ref. Q, i?c ~ Rq. As i?„ is tuned in 
such an experiment, one executes a typical trajectory A 
on the phase diagram of fig. 1 . The phase which exists in a 
real material like Gallium and is not captured by a pure 
model like model (1) is the Bose glass(BG) phase. So, 
instead of SC-BM-BI scenario along A, one would proba- 
bly have a SC-BM-BG(-Bf) scenario. Either way, in the 
BM phase, because of formula (25) the resistivity of this 
metal increases continuously from zero through a wide 
spectrum of values as i?„ is tuned through R^ as is seen 
in the experiments |qM| . We display the low temperature 
metallic resistivity from Ref. [||[fig.2 of this reference] as 
a function of i?„ in fig. 2 and observe the remarkable fit 
to the formula (25) with j/ ~ 2. This clearly points out 
that the low temperature metallic phase observed in Ref. 
fel was, in all probability, due to an incoherent motion of 



the Cooper pairs. The aforementioned value of the cor- 
relation exponent is consistent with the notion that Ga 
films were highly disordered samples |3^] and the Chayes' 
theorem [B^ along with it. 



V. HOW NEW IS THIS STATE OF MATTER? 

The question now arises whether the Bose metal dis- 
cussed in the previous sections is actually an adiabatic 
continuation of some known state of matter or if it is 
completely new. In this section, we argue that this state 
is actually a mathematical variant of a quantum disor- 
dered spin liquid with large spin and that it is a liquid of 
non fermionic variety. 

To show this, we first note that the algebra which con- 
trols the quantum mechanics of our problem is -E2, as 
discussed in section II. Now, E2 is a group contraction of 
S'0(3) ITs], i.e. as the radius R of the sphere on which 
50(3) operations are defined tends to infinity, iS'0(3) 
contracts to E2, viz. JxjR -^ —Py, Jy/R -^ Px, Jz — *■ L, 
where the operators L and P are as defined in eqn.(5). 
And, in the same limit, i.e. as the S0(3) angular momen- 
tum quantum number j = pR —^ 00, where p refers to the 
linear momentum quantum number of E2 , the irreducible 
representations of 50(3) map onto those of i?2 . Thus, the 
model (I) is a variation of a quantum spin model and the 
Bose metal obtained here, which is a quantum disordered 
phase, may be thought of as analogous to the frustrated 
spin liquids obtained within such a model in the limit of 
spin becoming large |p8| . 

A more intriguing question is the issue of how con- 
sistent the existence of a disordered phase as that of a 
Bose metal at T = is with the third law of thermody- 
namics. To this end, we calculate the low temperature 
specific heat of the metal under the same approximation 
as in section IV, i.e. non-interacting uncondensed bosons 
(vortices) moving in a transverse gauge field. The lon- 
gitudinal part of the gauge field is gapped because of 
the screening effects and hence do not contribute signifi- 
cantly to the specific heat p3|. However, the transverse 
gauge field representing the quantum fluctuations associ- 
ated with the transverse modes are gapless in this metal- 
lic phase |2^ (because there is no spontaneous symme- 
try breaking) and, hence, make a dominant contribution 
to the specific heat. The calculation for uncondensed 
bosons interacting with a transverse gauge field Go 33 
is analogous to the fermionic case for which an extensive 
literature exists |^^. Following Reizer |39|, we find at 
low temperatures (Appendix E) 



- AdTln{To/T) (with disorder) 



(26) 



This anomalous behaviour originates sheerly from the 
fact that the transverse gauge field undergoes dynam- 
ical screening, which leads to the presence of diffusive 
modes lu ^ ^iq'^, where n = 2(disordered) or 3(pure). 



The coefficient Ap scales as Ap 



2/3 



j: , Ad as Ad '^ uf, 
and To as To ~ ^/n/, where Uf refers to the free vortex 
density (Appendix E). This displays the non-fermi liquid 
behaviour of the Bose metal phase and verifies that there 
is no consistency problem with the third law. Clearly, the 
BM phase has more entropy than a normal fermi liquid 
at low temperatures. 

The foregoing temperature variation of the specific 
heat shows that there are weak singularities at T = in 
the metallic state. Thus, the BM state has the character 
of a critical point as T — > 0. This requires further study 
since the coefficients which enter the specific heat cal- 
culation would undergo strong renormalisation and the 
actual temperature dependence may be somewhat differ- 
ent. These arguments also suggest that there is no good 
length scale in the quantum liquid phase at T = 0. This 
can be seen by calculating the charge charge correlation 
function. Using equations (A4) (or, equivalently (B2)) 
and (E3), we obtained: 



< SUq.^Sn^q^^^ >? 



q 



-ia^ + bq^ 



(27) 



where the coefficients a and b may be read off from eqn. 
(E3). Thus, the excitation spectrum of the charge fluctu- 
ations is gapless and diffusive. As we pointed out in sec- 
tion IIB, because of the algebra of phase fluctuations, this 
does not contradict the fact that there is a length scale 
associated with superfluidity viz. the superfluid correla- 
tion length ^+ tied to the free vortex density nf. Thus, 
the BM phase is analogous to a gapless quantum spin 
liquid, very much like an RVB state. 



VI. DISCUSSION AND CONCLUSIONS 

In summary, we have argued that the superconductor- 
insulator transition in low— Tc superconductors should be 
viewed as a two order parameter problem and observed 
that it is described by two coupled XY models in (2-f-l)-d 
in the charge picture. This change from a single (2-|-l)-d 
XY model, as conventionally thought, is tied to the inclu- 
sion of an appropriate vortex mass term in the vortex pic- 
ture or non-local interactions in the charge picture. This 
leads to the interesting possibility of a novel disordered 
Bose metal phase distinct from the traditional supercon- 
ductor and Bose insulator phases. On the basis of our 
analysis of the model, we expect a superconductor - in- 
sulator transition when Von-sUe ^ Vr^.n. {n.n. — nearest 
neighbour), a case close to that of the Josephson Junction 
arrays (with B^xt = 0); and a superconductor - metal - 



insulator transition when Vo, 



K, 



a situation 



close to that of the granular superconductors (it may be 
helpful to note that Vi ^ (2e)^/Ci, Ci — capacitance) [Q. 
The physics of the problem seems to be controlled by a 
multicritical point. We find the properties of the Bose 
metal to be critical. It is inherently connected to the 



quantum spin liquid states and is strongly reminiscient 
of the gapless RVB phase. The resistivity of this novel 
metallic phase predicted by our calculations finds an ex- 
cellent match with the experiments on Gallium films [^ . 
Dissipation in the bosonic system at low temperatures, 
within our model, is hydrodynamical: it comes from the 
fact that a moving charge (vortex) dissipates as it moves 
through a background of charged particles (vortices and 
antivortices). The heat is carried away by the gapless 
transverse modes representing the plasmons. This source 
of dissipation translates into quantum fiuctuations in the 
charge picture. 

To our knowledge, so far there has been no metallic 
phase proposed in the phase diagram of the Bose locali- 
sation problem, within the scope of any physically realis- 
able model. This paper is a first attempt to propose this 
concept and establish some of the basic principles which 
underlie the existence of such a phase. As a result, many 
of the features which are tied to the current experiments 
have gone unaddressed. Some of these are: existence of 
disorder, long-range l/r interactions, non-zero magnetic 
field, finite temperature effects, etc. Also, in a realistic 
sample, an important source of heat bath in addition to 
the emission of transverse excitations is the electrons in 
the vortex core. This needs to be incorporated into the 
calculation. 

An important concern is whether the BM phase will 
survive in improved calculations, given that our calcu- 
lations have been mean field like. As we saw, the ex- 
periments argue in favour of the existence of this phase. 
Further, as we argued in section lid, the retardation ef- 
fects leading to the destruction of vortex superfiuidity 
and the binding of vortex - antivortex pairs, the features 
corresponding to the loss of charge order and the growth 
of phase order respectively, are controlled effectively by 
two separate parameters. Plus, as we found in section III, 
the constants which determine the curves LXM and ZYX 
in fig. 1 differ by about an order of magnitude because of 
the strong retardation effects felt by the vortices in this 
model. In view of these facts and the self consistent na- 
ture of the calculation, it is reasonable to think that the 
curve ZYX will lie above the curve LXM and intersect it 
at a multicritical point even if the fluctuations are taken 
into account, implying a separate BM phase E^. Either 
way, it will be very useful to check this via numerical sim- 
ulations. An important issue which confronts us at the 
moment is — what is the universality class of the crit- 
ical phenomena associated with model (1)? Although 
we have commented on this above, one needs to perform 
renormalisation group calculations or numerical simula- 
tions to answer this question completely. What is clear 
from our calculations, however, is that the critical be- 
haviour, both at T = and T = small but finite, implies 
a new universality class with a multicritical point which 
is very rich. 

Given the considerable amount of work done on the 
Bose-Hubbard models, one might wonder why this phase 
was not observed in the other theoretical constructs, par- 
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ticularly in simulations. In view of this, we would like to 
mention a few points, in addition to the idea of two order 
parameters, which distinguishes our work from the pre- 
vious ones. Firstly, we have worked in the large riQ— limit 
of the problem. Most of the current work on this model 
does not treat this limit consistently. Also, in this limit, 
the algebra which determines the quantum mechanics is 
qualitatively different from that of the hard core bosons, 
i.e. E2, and the chemical potential /x does not play a 
significant role (as found in section lib). In the small 
riQ— limit, fi plays a dominant role, especially in stabilis- 
ing the commensurate Bose insulator(CBI) phase (mostly 
called Mott insulator in the Bose literature). The weak- 
ening effect of /i as uq increases is clearly observable in 
the shrinking of CBI lobes with uq, evidenced by the 
perturbative and Quantum Monte Carlo calculations of 
Ref. |4^ . This provides additional support to the results 
obtained here in this approximation. Further, we also 
observed that the nearest neighbour interaction plays a 
crucial role in opening up the metallic phase in the sys- 
tem. There is no such phase with just on-site repulsion. 
This piece of physics has support from the RG flows con- 
structed by Fisher and Grinstein |1(| . 

On the supersolid phase: It is somewhat tempting to 
use eqn.(19) to the fullest extent and extend the curve 
ZYX beyond the point X along the dotted line XO. 
Eqn.(19) is really not good beyond the point X; but if 
we take it at face value, it seems to suggest a region 
OXM where both phase order and charge order exist, 
i.e. a supersolid phase. This has some support, if we 
look at the action (10). Here, as we noted earlier, when 
Vi <C Vb, 6*^ ~ 0, implying the presence of charge or- 
der in the system. This, however, corresponds to trivial 
charge ordering only, viz. 5ni ~ 0, i.e. Uj ~ no. This 
phase needs to be verified by better calculations. If this 
phase exists in improved calculations, the point X would 
be tetracritical. 

Before we close, we would like to make a short com- 
ment on the Mo-Ge system y], a situation where mag- 
netic field is present. As we argued in section lid, the 
physics of quantum phase fluctuations in 2-d supercon- 
ductors is a quantum mechanical generalisation of the 
KT scenario. We expect this to happen when magnetic 
field is present as well. Let us recollect the classical case 
for completeness. Here, dislocation - antidislocation pairs 
are induced in the vortex lattice which unbind as a result 
of thermal fluctuations ^^ . At low temperatures, quan- 
tum effects dominate and we expect the unbinding result- 
ing from quantum fluctuations instead. But, as pointed 
out in this paper, their kinetic energy should receive two 
contributions — from the zero point motion and the ac- 
tion of a heat bath. The metallic state is realised when 
the dislocations in the vortex lattice(more precisely, vor- 
tex glass) phase move under the action of a heat bath. 
This is very much along the lines of what had been sug- 
gested in Ref. |Q. Following our discussion in section V, 
we suggest specific heat measurements be made to spot 
any non-fermionic behaviour of this metal-like phase. We 



plan to take up some of the unresolved issues mentioned 
in this section in a future publication. 
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APPENDIX A: THE DUALITY 
TRANSFORMATION 

We start from action (9), but with the transformation 
rrii — > e~*'^'"*mi undone, so that 

S = iY^ m,{V^(j),) + VoAt J2 "1? 

i i 

+ ViAr ^(mi -f TOi+a)2 - JAr ^ cos(Va0i)- (Al) 

ia ia 

We now use Villain transformation on the J term. For 
small values of J At, one has |34] 

JArV cos(Vc0,) -'"(T^)y^ "L + 'V 0.(VQnic) 

Integrating out the phase degrees of freedom cpi , one ob- 
tains — 

{nii} {Tii„} 



ViAT^{mi + mi+af 



supplemented with the constraint 

Vrmj +Vania == 0, 



(A2) 



(A3) 



(with the summation over a implied). Defining n^ — 
{mi, nf ),(// = 0, 1, 2), one obtains the constraint (A3) as 



V^< = 0. 



This allows us to define an integer gauge field Af , via 



— £ fii/p V f-^i • 



(A4) 



A purist might object to our not using the shifted lat- 
tice operators A^_~ in eqn.(A4) (l3|. This does not 
matter for our case since the model involves a single 
gauge field. Now, transforming over to the Fourier space 
{rrii = ^7= J2q,uj e''"^""!,,^, with q= {uj,q) and Xi refer- 
ring to the i*'' site on the space-time lattice, with Nq the 
total no. of lattice sites), one obtains. 
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with 



S^ln{-—)Y,[eAl^A\_^ + G-\uj,q)Al^A%^_^ 



(A5) 



with G ^{u!,q) given by eqn.(12) and A^ ^ refer to the 
fourier transform of A'^. Here q and a) refer to lat- 
tice momentum and frequency respectively, and q^ — 
Eo.Kaiq)K^iq), with ;^„(<7) = i(^_ e"'?"), i^„(g) = 
l{e^i^^l)^. And, (q-Q)' ^j:^K^{q-Q)K^{q-Q) 
in eqn.(12). Eqn.(14) is complete with the gauge condi- 
tion, 



V„^" = 0. 



(A6) 



Now, we introduce the vortex variables jf via the Poisson 
summation formula llSl — 



bTl 



(A7) 



Integrating out the gauge field degree of freedom, one 
obtains eqn.(ll). 



APPENDIX B: CONTINUUM VERSION 

The features of the vortices discussed in section lid are 
not special to the lattice case, but exists in the continuum 
case as well. This gives us a more general frame in which 
to study the model. The considerations in the continuum 
case start from the action p3,M — 

S^ I dTd'x[{m/2p)\J\^ +\ [ d?yV{x~y)5p{x)5p{v) 



2T:ij^A^ 



(Bl) 



where J refers to the charge current and Sp the charge 
fluctuations, with the gauge field A^ given by — 



J'^(x) = (V,J) 



' llUpCfu-ri-p 



(B2) 



Also, jp — (py, Jy) refers to the vortex three current den- 
sity, m — mass of the Cooper pairs, and p — average 
Cooper pair density. We have not written down a Mag- 
nus force term on vortices due to a background charge 
condensate, because this is not contained in the lattice 
model (1) discussed in the paper. However, this may be 
important in a uniform system like Mo-Ge H]. Work- 
ing in the transverse gauge V.A = 0, we obtain after 
integrating out the gauge field — 



1 .n 






(B3) 



where the Green's function G(^lj, k) is given by — 

G"i(tj, fc) = w^ + {p/m)k^Vk (B4) 

where Vk is the Fourier transform of the potential V{r). 
From here onwards we consider the screened Coulomb 
potential with 



Vk = Voa/VfcM^ 

In the limit tj,fc — > and to <^ Csk (with Cg being the 
plasmon velocity), the Green's function splits up into 
a singular part Gs(td,k) and Go just like in the lattice 
model — 

G{uj,k)~Gs(uj,k) + Go, (B5) 

where 

Gs{iu, k) = 1/(0.2 + elk"), Go = l/{2pVoaym), 

with 

Cg = pVo/m. 

Thus, from (B3) and (B5), we obtain — 



S = 2tt 



2 P. 

m 



7 .Jk, 



-7° 



^^'^Y.~^kJGsiu;,k) + G„)~f_ 



k,—uj 



(B6) 



Thus, comparing eqns.(B6) and (14), we see that the lat- 
tice and continuum actions of phase fluctuations, within 
the scope of our model, have identical low frequency and 
long wavelength limits. 

Correspondence between the lattice version and the 
continuum version is established via the following iden- 
tifications — 



m/ p 
pVo/m 
Go - ll{2pVoa'lm) • 



2ln{2/JAT) 



(B7) 



a? 



cj = {Vq + 8Vi)AT/ln{2/JAT) (B8) 
Go = b^/8{b^ + c^f (B9) 



Following Ref. |4^, one can easily show that eqn.(B3) is 
equivalent to the following Hamiltonian for the vortices 



Hy — 



1 



2m, 



'S^ipi + 2TTqiA{xi)Y + 2'k'^— V" q.,qAn \ n ~ r, \ 
^ — ' m -"^ — ' 

(BIO) 



where the summation is over all vortices and antivortices 
and the gauge field A has the spectrum determined by 
Gj^{uj,k) as in eqn.(15). Here qi = ±1 refers to the 
charges on the vortices. Here my refers to the vortex 
mass. 
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APPENDIX C: AN ESTIMATE OF G(0) 

From section Ilia, we have — 
1 r duj r (Pq 



=<°'^^/.j/: 



1 



(Cl) 



(27r)2 (1 - cosuj) + A2 
where 

A^ — \2— {cosQx + cosqy)][KQ + 2ti\{2 + cosqx + cosqy}] 

(C2) 

with 

Kq = (? and k\ = 6^/8, 

where 6^ and c^ have been introduced in section lid. 
From (Cl), we obtain 



=(°'^^£(^^<"' 



with 



F{<1) 



VM2(2 + A2) 



(C3) 



(C4) 



Important contributions to (C3) come from the low mo- 
mentum region. Thus, 



F{q) ~ 1/^2(^2 + 8^2)^2 - {cosq^ + cosqy) (C5) 

Combining cqns.(C3) and (C5), we are led to the condi- 
tion (17), with 



Jo 



^/^ <P, 



\J sin^x + siv?y 



(C6) 



APPENDIX D: SELF CONSISTENT 
FUNCTIONAL APPROACH 

In this appendix, we give a brief review of the self- 
consistent functional approach Gsl discussed by loffe et. 
al. 12^. This is essentially a hydrodynamical kind of 
approach. We start with the following electromagnetic 
Lagrangian of bosons in imaginary time (in 2 — d) — 



i^B = i^B + ^B + i^B- 



(Dl) 



Here 



i°B = ^* ^ + T5^ I W P + ^P(r)/n I r - / I p(r') 
OT 2m 2 



Mi^r 



{D2a) 



Lb = -^a.[^pvr-r^i^] + 2^a' 1^1' - (^26) 



L% = ^iidraf + c2(V X df] - {D2c) 

where ip,ip* refer to the Bose fields and p — tp*ip. We 
work in the transverse gauge V.a = 0. The dimension- 
less coupling constants determining the strength of the 
Coulomb repulsion and transverse gauge field, viz. ac 
and ag respectively, are defined as — 






(D3) 



where n is the average density of bosons. Comparison of 
eqn.(D2) with eqn.(15) allows us to identify — 



5^- 



2 2 

C ^^ C, 



2Tr^/ln{2/JAT) 



and, 



m/2n ^ T:^Go/ln{2/JAT) 



(D4) 



The last identification follows from the fact that the 
weight of the zero point motion term | j p is m/2n. 
Using eqns. (D3) and (D4), we are led to eqns.(18) Q. 

The calculation of the depletion of the Bose conden- 
sate from the action (D2) proceeds in two steps. The 
first involves determining the effective functional of the 
superfluid Bose system in the absence of a transverse 
gauge field. And, the second involves determining the 
effect of the gauge field on this system, i.e. calculating 
the change in the superfluid density. 

The effective functional of the superfluid Bose system 
(with a = 0) is obtained by expanding the Bose field in 
terms of the slow and fast modes, and integrating out 
the fast modes i4ql . Galilean invariance of this system 
dictates Ug — n ||26|] . The leading terms in the effective 
action are (for small values of a) — 



So 



d xdrl 



'2m 



Vcl) 



indrCJ) - 






1 



8mn 
1 



Vtt 



(D5) 



where </> and p{x, t) — n + n{x, t) are the phase and the 
amplitude of the slow mode -00 — \/p{x,T)e^'^^^-''^\ We 
now work in the real time and integrate out the phase 
variables. This leads to an action solely in terms of the 
density fluctuation variables tt — 



■ J d'^xd^x'Tr(x)C^'^{x-x')7T(x') 



where 



C(^,g) 



ru^ jm 
up- — ((p I2wi)'^ — ng^ jm + ib 



(D6) 



and iPx = Pxdt. 
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To calculate the change of the superfluid density Ug = 
n due to the action of the transverse gauge field, we focus 
on the diamagnetic term in (D2b) 

'^-'ro^o^-^An + ^y (D7) 

The weight of the a^ term determines the superfluid den- 
sity. One now integrates out the 7r-degrees of freedom. 
The coupling term in (D7) leads to an cP — cP interaction 
term. One then (a) decouples this term in a mean field 
approximation and (b) rewrites the gauge field in terms 
of fast and slow modes , ai and ap, and integrates out 
the fast modes. The weight of the Qq term then gives the 
self-consistent equation for the superfluid density rig, 



ris = n-\ / (fxC{x)D{x) 

m J 

where D{x) is the gauge field correlator ■ 

■,2 



DiiO.q) 



9 



llP — [(?<f- + UsQ^ jni) -\- iS 



(D8) 



(D9) 



Eqn.(D8) leads to the phase diagram in a — ag plane, 
mentioned in section Illb. 



APPENDIX E: SPECIFIC HEAT CALCULATION 

To calculate the specific heat, we consider the La- 
grangian (D2) without the log interation, which is 
screened in the metallic phase. Connection with the orig- 
inal model is established via the mapping (D4) JS^I. We 
first rescale the transverse gauge field a ^ ga so that the 
effective charge g appears in the eqn.(D2b) rather than 
(D2c). Thus, the unperturbed gauge propagator is 



Dlfjiq^UJn) ^Kafiiq) 



1 



UJt, 



c?q'^ 



(El) 



where w„ = 2TmT{n =integer) are Matsubara frequen- 
cies and 

Kafj{q) = Oaf3 ^. 

To obtain the contribution of specific heat from the trans- 
verse modes, one integrates out the particle degrees of 
freedom in a one- loop approximation p9[ . This leads to 
the retarded polarization function — 



ni?(<7,^)=-ifa/3(g)5'[- 



2n n 

I 

m vsq 



+ c^oq^] (E2) 



where xd = ns(— /i)/247rmc^ and mvg/2 —\ ^b |, with 
nsiO = l/(e'^^ - 1). The structure of eqn.(E2) is quite 
independent of the statistics because the integrals which 
enter the calculation are of the form J_ d£,n{£,) and 

JT d^W- "^^^ renormahsed gauge propagator is then 
given by Dyson's eqn. — 



D = Do + DqUD 

In the quasistatic approximation fl ^ cq, we have the 
retarded gauge propagator as 



tiR) 



D'^'^{q,^)^K^p{q)[l/{-z 



2ng^ n 
m vgq 



+ cV)], (E3) 



where c^ = c'^(l+g'^XD)- Integrating out the gauge fields 
leads to the free energy — 



PF = --TrlnD{q, iVn) 



(E4) 



The specific heat is obtained from C — TdS/dT. The 
details of the calculation may be found in Rcf . |^ . For 
the disordered case, the Green's function is written as 
G{p,LUn) — \/{iuJn — ^p + i /2Tsgn(ujn)) , with T the scat- 
tering time for the vortices from the impurities (which 
appear as a static gauge field in the vortex picture) P4| . 
In this case, the factor X/vsq appearing in the frequency 
term in (E2) and (E3) gets replaced by 2t. The coeffi- 
cients Ap, Ad and To in eqn. (26) are approximately given 
by- 



Ar, 



{2ng^lmvBC^f'^l-K 



\ r^, 2/2-2 

Ad ~ ng T/TT mc 
To ~ qc^mc^ /4:ng'^T 



(E5) 



where qc denotes an upper momentum cutoff ~ 1/Coj 
with ^0 the Cooper pair size. 



APPENDIX F: COMMUTATION RULE (6C) 

In this appendix, we discuss how the commutation rule 
(6c) comes about. From commutators (6a) and (6b), we 
observe that P and P^ are ladder operators. Hence, in 
the angular momentum representation (L | flm >= m \ 
f},„ >) 0,|l|, 

oo oo 

p^i Y^ I f7„_i >< a„, |« i ^ I f^m-i >< n„, I 
P^ ^-i Yl I ^"+1 >< ^™ i~ ~*X1 1 ^™+i >< ^™ I 

-{no} -oo 

where {no} = uq for integer no, and the floor or ceiling 
of no (appropriately taken) otherwise. The approximate 
expressions written on the right hold only when no S> 1. 
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FIG. 1. A schematic phase diagram for model (1). The var- 
ious phases are demarcated by solid boundaries. X denotes 
the multicritical point. The dashed lines (OA, OB, OC) de- 
note typical trajectories executed as the normal resistance R„ 
of a superconducting film is tuned. Please see the text for de- 
tails. 



FIG. 2. The resistivity of the metallic phase as a function of 
R„ in Ga film experiment [5] (the diamonds are the datapoints 
from fig. 2 of this reference) and the best fit of eqn.(25) (the 
dashed curve). (We have not displayed those points of fig. 2 
in Ref.5 which are in the insulating regime and, expectedly, 
deviate from the best fit shown here.) 
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